Which Kubo formula gives the exact conductance of a mesoscopic disordered system? 
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In both research and textbook Uterature one often finds two "difi'erent" Kubo formulas for the zero- 
temperature conductance of a non-interacting Fermi system. They contain a trace of the product of 
velocity operators and single-particle (retarded and advanced) Green operators: Tr{vx & Vx G") or 
Tr{i)x ImG Vx^ra G). The study investigates the relationship between these expressions, as well as the 
requirements of current conservation, through exact evaluation of such quantum-mechanical traces 
for a nanoscale (containing 1000 atoms) mesoscopic disordered conductor. The traces are computed 
in the semiclassical regime (where disorder is weak) and, more importantly, in the nonperturbative 
transport regime (including the region around localization-delocalization transition) where concept 
of mean free path ceases to exist. Since quantum interference effects for such strong disorder are 
not amenable to diagrammatic or nonlinear a-model techniques, the evolution of different Green 
function terms with disorder strength provides novel insight into the development of an Anderson 
localized phase. 
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in practice, as w« 
in both research! 



At first sight, the title of this paper might sound pei« 
plexing. What else can be said about Kubo formulatl 
after almost a half of a (last) century of explorations 
1 as through nyyaerous re-derivations 
and textbookoQ literature? Kubo 
linear response theory (KLRT) represents the first full 
quantum-mechanical transport formalism. It connects 
irreversible processes in nonequilibrium to the thermal 
fluctuations in equilibrium [fluctuation-dissipation theo- 
rem (FDT)]. Therefore, the study of transport is limited 
to the nonequilibrium states close to equilibrium. Nev- 
ertheless, the computation of linear kinetic cocfflcients is 
greatly facilitated since final expressions deal with equi- 
librium expectation values of relevant physical quantities 
(which are rrmch simpler thai|Lthe corresponding nonequi- 
librium onesO). It originatedQ from the Einstein relation 
for the diffusion constant and mobility of a particle per- 
forming a random walk. i-. 

Until the scaling theory of localization,!] and ensu- 
ing computation of the lowest order quantum correc- 
tion, weak localizationQ (WL), to the Drudc conductiv- 
ity, it almost appeared that microscopic and complicated 
Kubo formulation of quantum transport merely served 
to justify the intuitive Bloch-Boltzmann semiclassical ap- 
proachB to transport in weakly disordered {kpi ^ 1, kp 
is the Fermi wave vector and £ is the mean free path) con=j 
ductors. Furthermore, the advent of mesoscopic physicsEl 
has led to reexamination of major transport ideas — in 
particular, we learned how to apply properly KLRT to 
flnite-|aize systems. Thus, the equivalence was estab- 
lisheda between the rigorous KubOpEprmalism and heuris- 
tically founded Landauer-Biittikert3 scattering approach 
to linear response transport of non-interacting quasipar- 
ticles.B This has emerged as an important tool in for 
studying mesoscopic transport phenomena, where system 
size and interfaces through which electron can enter or 
leave the conductoij play an essential role in determining 
the conductance.BEj 



This study presents an exact evaluation of two differ- 
ent Kubo-type expressions for the linear conductance of 
a mesoscopic disordered conductor. Both expressions are 
frequently encountered in research as well as textbook lit- 
erature, and are displayed below. They consist of a trace 
(or linear combination of such traces) over the product of 
velocity operators Vx with retarded and advanced single- 
particle Green operators G""'" = [E — H ± zO'*']""'^, like 



Tr 



where H is an equilibrium Hamil- 

tonian (in the spirit of FDT, it contains random and 
confining potentials, but not the external electric field), 
and velocity operator is defined by ihv = [r,H]. These 
quantum-mechanical traces are computed here. Figs. ^ 
and ^, in the site representation (i.e., using real-space 
Green functions) defined by a lattice model, such as tight- 
binding Hamiltonian 



H = ^em|m>(m| 



E 

(m,n) 



tmn|m)(n| 



(1) 



on a hypercubic lattice N"^ of size L — Na {a being the 
lattice constant). Here tmn is the nearest-neighbor hop- 
ping integral between s-orbitals (r|m) = tp{r — m) on 
adjacent atoms located at sites m of the lattice (tmn = 1 
inside the sample defines the unit of energy). The dis- 
order is simulated by taking random on-site potential 
such that Em is uniformly distributed over the interval 
[—W/2,W/2], which is the so-called Anderson model of 
localization. I emphasize the requirements of current con- 
servation throughout this analysis, which will allow us to 
understand the features of different trace expressions in- 
troduced above. |-||— ■ 

The mesoscopic naathods (mesoscopic Kubca-llil or, 
equivalent, Landauerlij formula) make it possible to get 
the exact zero-temperature (i.e., quantum) conductance 
of a flnite-sizc sample attached to semi-infinite disorder- 
free leads. Although KLRT is a standard formalism for 
introducing the many-body physics into the computation 
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FIG. 1. Different terms in the Kubo formula for the 
two-probe quantum conductance of a single finite-size sam- 
ple modeled on a simple cubic lattice 10^ by an Anderson 
model with disorder strength W = 2 [upper panel — single 
sample in the semiclassical transport regime] or = 7 [lower 
panel — disorder- averaged over 50 samples in the nonpertur- 
bative transport regime kpi ^ 1]. The full Kubo conductance 
(thick solid line) is given by the sum of terms defined in Eq. (^) 
[thin solid line] and Eq. (m [dashed line], G = Gra + Grr- The 
respective traces in these expressions are performed only over 
the states residing on the first two planes inside the sample. 
The dotted line in the upper panel represents G^ra'^ obtained 
by tracing over the two planes deep inside the left lead (at 
the distance 10a away from the sample). 

of transport coefficientsj^ here the focus is on the 
transport properties determined by scattering of non- 
interacting (quasi)electrons on impurities. The "old" 
Kubo formulaEj for the macroscopic volume-averaged 
longitudinal DC conductivity at zero temperature {E = 
Ep in all formulas below, Ep being the Fermi energy) of 
a non-interacting Fermi gas described by a single-particle 
Hamiltonian H is given by 



(Txx — — — Tr 



VxS{E-H)vx6{E-H) 



(2) 



FIG. 2. Different terms in the Kubo formula G = Gra+Grr 
for the two-probe quantum conductance of the same finite-size 
conductors as in Fig. ^, but with respective traces in these 
expressions performed over the site states inside the whole 
disordered sample. 

j = /drj(r)/ri to the spatially-averaged electric field 

oo 
im- 



fiE, where thermodynamic limit ft 



where factor of two accounts for the spin degeneracy. The 
Kubo conductivity relates the spatially averaged current 



(while keeping the impurity concentration finite! Jis 
plied to get the unambiguous intensive quantityll3'lj (and 
well-defined steady state). For electrons in a random po- 
tential further averaging should be pcrfbimcd over the 
disorder to get cr as a material constant .EZl On the other 
hand, quantum corrections to the conductivity are non- 
locaB on the scale of the dephasing length 3> i- This 
invalidates the concept of local quantities, like conduc- 
tivity, in mesoscopic samples, which are smaller than 
and thereby effectively at T = 0. Therefore, mesoscopic 
transport has to be described in terms of sample-specific 
quantities, like conductance, which describe a given sam- 
ple measured in a given mannero (i.e., more generally, 
conductance coefficients Ip = dpqVq in the Ohm's law 
for a multi-probe geometry, where several leads are at- 
tached to the sample to feed the current Ip or measure 
the voltages Vq), or alternatively, non-local conductivity 
tensor introduced below. Switching to conductance leads 
to a following Kubo expression 



^ 4e2 1 
1 



hvx Im G hvx Im G 



ImG = -(G"' - G") = -ttS(E - H). 
2i 



(3a) 
(3b) 



Here the definition of retarded (r) or advanced (a) single- 
particle Green operator G""'° — [E ~ H ± zO+]~^ re- 
quires a numerical trick to handle the small imaginary 
part iO^, which then spoils the prosDeC|l_pf obtaining the 
exact zero-temperature conductance J130 Once the semi- 
infinite clean leads are attached to the finite sample (at 
planes 1 and N along x-axis for a two-probe geometry, 
Fig. ^), the "self-energy" E'"'" = E^'" + S^", arising from 
the "interaction" with the leads (L-left, i?-right), pro- 
vides a well-definedr-imaginary part in the definition of 
the Green operatorsEJ 



(4) 



The Green function G'"^''(n, m) = (n|G'''°|m) describes 
the propagation of electron between two sites inside an 
open conductor [L^ = L vn the two- probe geometry). 



The self-energy terms are givenEJ^EJ by E£ ^ (n, m) = 
(^c/2t£)5£_^(ns,ms) with 5£_^(ns, ms) being the sur- 
face Green function of the bare semi-infinite lead between 
the sites ns and 1115 located on the end atomic layer of 
the lead (and adjacent to the corresponding sites n and 
m inside the conductor). It has an imaginary part only 
for \E\ < G^L, which means that G{Ep) goes to zero at 
band edge of a clean lead \Ep\ = G^l because there are 
no states in the leads beyond this energy which can carry 
the current. Here the leads will be described by the TBH 
with £m = and tmn = tv, the hopping between the sites 
in the lead and the sample is tmn = tc, ^ illustrated in 

I use the term "mesoscopic Kubo formula" for Eq. (|^) 
with the Green operators r-ZQ) plugged in. This for- 
mula ia_exactly equivalentEl to a two-probe Landauer 
formulall3 for the conductance measured between two 
points deep inside macroscooic reservoirs to which the 
leads are attached at infinity.Ej Thus, it is conceptually 
different from the "plain" Kubo formula (||) which fol- 
lows from combining the conductance of smaller parts 
G = aL'^^^/L, thereby implying local description of 
transport which becomes applicable only at sufficiently 
high temperatures. Such mesoscopic formulas provide 
means to compute the quantum conductance — a sample 
specific quantity which takes into account the finite sys- 
tem size, measuring geometry, arrangement of impuri- 
ties, non-local features of quantum transport, and can 
describe ballistic transport (where local relation j = crE 
does not hold). Attempts to use the original Kubo 
formulas, xin finite samples, throughout premesoscopic 
historylia^Ej of the Anderson localization theory, were 
thwarted with ambiguities, which can be traced back to 
the general questions on the origin of dissipation.^ This 
stems from the fact that stationary regime cannot be 
reached unless the system is infinite or coupled to a 
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FIG. 3. A two-dimensional version of the actual 3D model 
of a two-probe measuring geometry employed here. Each site 
hosts a single s-orbital which hops to six (or fewer for surface 
atoms) nearest neighbors. The hopping matrix element is t 
(within the sample), th (within the leads), and tc (coupling 
of the sample to the leads) . The leads are semi-infinite and 
are considered to be connected at ±cxd to "reservoirs" biased 
by the potential difference eV = /il — I-Ir- 

thermostat. From the practical point of view, leads 
make the system infinite by opening the sample, and 
therefore eliminating the technical obstacles ia handling 
of the discrete spectrum of finite samples.c3 Further- 
more, the usage of semi-infinite leads allows us to by- 
pass explicit modeling of reservoirs in the computation 
of conductance since "hot" electrons which escape into 
the leads (due to the broadening of energy levels) do not 
come back in a phase-coherent fashion. The concept of 
reservoirs wasi-sdways essential part of Landauer 's sub- 
tle arguments. r4 They provide dissipation and therefore 
the steady statc.EZI However, the computation of conduc- 
tance, as a measure the dissipation, involves only a con- 
servative Hamiltonian of noninteracting electrons scat- 
tered on impurities (i.e., when computing linear conduc- 
tance of a mesoscopic system one usually does not deal 
explicitly fmth electron-electron and electron-phonon in- 
teractionstiJ) . 

The conductance computed from the mesoscopic quan- 
tum expression is exact, but characterizes the whole sam- 
ple+leads system in the spirit of quantum measurement 
theory since leads can also be papLidered as the "macro- 
scopic measuring apparatus" However, for disor- 
dered enough sample (and not too narrow leads or too 
small tc) the conductaii(je-|is determined mostly by the 
disordered region itself.E§Eil This exactness makes it pos- 
sible to compute the transport properties in both semi- 
classical regime (where Boltzmann theory and perturba- 
tive quantum corrections are applicable since kpi ^ 
as well as in the nonperturbative transport regimeE3 
where semiclassical concepts, like £, loose their mean- 
ing. Although the distinction between kp£ ^ 1 (where 
semiclassical theory, including the perturbative quantum 
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corrections, breaks down) and the criterion G ^ 2e^ /h 
for the locahzation-delocahzation (LD) critical point has 
been known since the scaling theory of localization,!] it is 
not uncommon practice to find these two different bound- 
aries of transport regimes confused. As the disorder is 
increased, a sample goes from the semiclassical transport 
regime, through a vast region between kpt ~ 1 and LD 
transition Ifi.g., 6 < < 16.5 at Ep = in the Ander- 
son modelo), and finally enters into a localized phase. 
Since the nonperturbative transport regime lacks small 
parameter required by present analytical schemes, the 
numerical techniques employed here are the only way to 
gain insight into the quantum effects beyond the lowest 
order corrections, like WL, or resummation of all such 
perturbative quaiatjHH correction within the nonlinear a- 
model formalism E3^Eil 

Anflfegfi exp ression is often encountered in the liter- 



for both conductance of finite-size sys- 



aturi 

terns and conductivity of infinite systems. It gives 
the conductance through an apparently different trace. 



G oc Tr 



hVr G'' hVr G" 



While some textbooks quote 

this only as a conveni|ept||a,pproximation to the disorder- 
averaged full formula,alljia sometimes the claim goes fur- 
ther to say that such trace is equivalent to the Landauer 
formula, and moreover it can be evaluated at any cross 
section inside the disordered region, thus relaxing the re- 
quirement of perfect le^ds "as the weakest point of the 
Landauer formalism" E3 Namely, G' G"' or G°'G'^ terms 
can be reduced to a single Green function^ia a Ward 
identity in weakly disordered conductors Ej'u and are 
therefore related to the density of states. They are "aban- 
doned" in the limit kpi ^ 1 since they do not generate 
interesting contributions to WL or mesoscopic fluctua- 
tions.liB In fact, both the diffusion modes of the Diffuson- 
Cooperon diagrammatic perturbation theory and differ- 
ent •vtersinns of the nonlinear cr-model (NLSM) are de- 
rivedc^Eil by considering only the term containing the. 
product G^'G". The NLSM is a quantum field theoryO 
of weakly disordered mesoscopic conductors where diffu- 
sion modes of the diagrammatic perturbation theory play 
the role of soft modes responsible for long-range spatial 
correlations of local density of states, mesoscopic fluctu- 
ations of glohal quantities, and nonlocal corrections to 
conductivity.^ It makes it possible to handle the break- 
down of pcctiirbation theory due to the proliferation of 
such modesLj by, summarizing all WL-type corrections 
to conductivity.E3 This then justifies the phenomenologi- 
cal one-parameter scaling theoryQ (if not rigorously, then 
at least qualitatively), and explains the LD transition in 
2 -\- e dimensions where Anderson localization occurs at 
weak disorder kpi:s> 1. 

To remove possible confusionjHj it should be empha- 
sized that conductance coefficients gpq, obtained by inte- 
grating Kubo non-local conductivity tensor g(r, r') over 
the cross sections in the leads p, q 



dSp ■ (T(r,r') • dS 



9' 



(5) 



also contain Im G for p ~ q, whilC|-|for p =/= q only the 
terms involving G^G"" are non-zero.l3 The cross sections 
Sp and Sq are to be chosen far enough from the sample, 
where all evanescent modes have died out. This not only 
provides the rigorous foundation for the Landauer for- 
malism, but also clarifies some subtle, points in the Kubo 
formalism (like disorder averaging^ and independence 
of linear transport ppoperties from the nonequilibrium 
charge redistributionEZi) . By writing the full Kubo for- 
mula (^) as a sum G = G, a + Grr , I introduce the "Kubo 
conductance terms" 



r _ 1 



2e^ 1 



hv^ G"" hv^ G" 
hvx G^ hvx G^ 



(6) 
(7) 



Obviously, if in some transport regime conductance can 
be obtained from the trace in G^a, the other term G^r 
has to vanish, at least approximately. 

Before embarking on the direct evaluation of these 
expressions for a conductor described by TBH (|^), the 
crucial point is to understand the way (i.e., space of 
states |m}) in which the traces should be performed, 
Tr(. . .) — J2mi^\i- ■ Ol"^)- Naively, in the site represen- 
tation it would appear that trace in formula (^ should 
include site states inside the whole sample. However, 
once the current conservation V • j(r) = is invoked this 
becomes extraneous. All Kubo conductivity or conduc- 
tance formulas stem from the more fundamental quantity 
in KLRT, the non-local conductivity tensor relating local 
current density to the local electric field, 



j(r)= y'dr'g(r,r')-E(r') 



(8) 



This tensor is obtained as a response to an external field 
only since corrections to the current due to the field of-jp-i 
duced charges go beyond the linear transport regimeEj'E^I 
(i.e., one does not have to engage in a much more cum- 
bersome task of finduag the response to a full electric field 
inside the conductoro) . In application of KLRT to meso- 
scopic systems, ^(r, r') is a sample-specific quantity, i.e., 
defined for each impurity configuration and arrangement 
of the leaddp, and is in fact non-local even after disorder 
averagingO Because it is not directly measurable, some 
volume averaging is needed to get the quantities that can 
be related to experiments 



G=-^y" drdr'E(r)-g(r,r')-E(r'). 



(9) 



Here V is the bias voltage {V ^ in the linear transport 
regime), e.g., in the case of two-probe geometry eV = 
Mi ~ M-Rj where leads are in equilibrium with two 
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FIG. 4. Different terms in the Kubo formula G = Gra+Grr 
for the two-probe quantum conductance of the finite-size 
conductor modeled on a lattice 10^ with diagonal disorder 
= 15 (whole band becomes localized at Wc — 16.5). 
Respective traces in these expressions are performed over 
the site states inside the whole disordered sample and dis- 
order-averaged over 50 realizations. 

macroscopic reservoirs characterized by a constant 
chemical potentials and /i^ (Fig. ^). Although this ex- 
pression contains the local electric field E(r) = ViJ.{r)/e 
inside the conductor, because of current conservation en- 
tailing V • g(r,r') — g(r,r') • V' = (which is a special 
case, in ihe absence of magnetic field, of a more general 
theoremP) , the conductance will not depend on this field 
factorsEJ. In the case of TBH (^ with nearest-neighbor 
hopping the expectation value of the velocity operator in 
the site representation (m|{}j:|n) — {i / K) t-^i-a^ {nix — n^) 
is non-zero only between the states residing on adjacent 
planes (technical details of a route from Eq. (|^) to the 
trace involving velocity operator are meticulously cov- 
ered in Ref. Thus, the minimal space choice here is 
defined by taking E(r) to be non-zero on two adjacent 
planes, while thc-standard textbook assumption of a ho- 
mogeneous fieldHila |E| = V/L throughout the sample 
leads to Eqs. (||), (|3[). Blind tracing over the whole sam- 
ple would give simply the conductance multiplied by the 
square of the number of pairs of adjacent planes, mean- 
ing that such trace should be divided by {N — l)^a^ to 
get G [when trace is performed only over the arbitrary 
two adjacent planes, in Eq. (||) is replaced by a?]. 
The physical content of this statement is simple: current 
/ = GV is the same on each cross section. Therefore, 
the conductance depends only on the total voltage drop 
over the sample, and not on the local current density and 
electric-field distributions. 

At this point, it is worth mentioning that another ex- 
pression is frequently employed in the real-space compu- 
tational practice. It stems from the linear response limit 



of a formu 
interacting 



-derived by the Keldysk-technique (for non- 
3 or interacting systemsCj) 
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FIG. 5. Different terms in the Kubo formula G = G,.a+Gr,- 
for the two-probe quantum conductance of a clean finite-size 
conductor modeled on a lattice 3^ (em = Q, t = t-^ — tc ~ 1) ■ 
The conductance as a function of the Fermi energy changes in 
steps corresponding to the number of open conducting chan- 
nels (defined by 9 quantized transverse propagating modes). 



G 



4e^ 



Tr ImSLG^TvImSflG^i 



(10) 



Here G\j^, are the submatrices of the full Green 

function G"'''^(m, n) whose elements connect the layer 1 
to layer N of the sample. Therefore, only a block N'^ x N'^ 
of the whole iV'^ x matrix G"'(n, m) is needed to com- 
pute the conductance. Although such "partial" knowl- 
edge of the whole Green function and the trace over ma- 
trices of size N'^ in Eq. ( [l0| ) are different from the cor- 
responding counterparts required in application of the 
mesoscopic Kubo formula (where minimal trace goes over 
2N'^ X 2N'^ matrices), the final result for the conduc- 
tance is the same. Positive definiteness of the operators 
—2lmT,Lji makes it possible to find their square root 
and recast the expression under the trace of Eq. ( p^ as 
a Hermitian operator. The expression then looks like 
the two-probe Landauer formula involving a transmission 
matrix t 



G=-Tr(ttt) = ::^ 



t = 2V-Im SlG^ 



N 



-ImS; 



(11) 



(12) 



or transmission eigenvalues r„ when the trace is evalu- 
ated in a basis which diagonalizes tt^. Moreover, it is 
more efficient from a practical point of view since Green 
function techniques used to evaluate Eq. (|l^) do not re- 
quire to know exact asymptotic eigenstates in the leads 
(as is the case with "original" Landauer formulation re- 
lying on the knowledge of the scattering basis of wave 
functions defined within the asymptotic regions of the 
leads). This becomes sine qua non when computing the 
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conductance of systems with complicated "Leads", e.g., 
like in the case of atomic-size point contacts.EJ 

The trace expressions like and (|^) do not conserve 
the current. Therefore, the suggested strategy is to evalu- 
ate Gra and Grr by tracing over the states located on two 
planes inside the sample (or inside the leads), as well as 
over the whole sample (and see how close Gra can get to 
G) . Different types of conductors are chosen for this eval- 
uation: weakly disordered with W = 2 (e.g., ^ w 9a at 
Ep = 0), and strongly disordered conductors with W = 7 
(for which unwartanted use of the Drudc-Boltzmann for- 
mula would givec3 £ < a). The exact computation of 
these traces is shown in the upper panels of Figs. |^ and 
H for a single impurity configuration oiW — 2 (here Gra 
is also computed by tracing over two planes deep inside 
the leads, which corresponds to integrating a{r,r') over 
such cross section, as discussed above). The lower pan- 
els of Figs. 1^, H plot disorder-averaged quantities over an 
ensemble of impurity configuration for disorder strength 
W — 7. In both ways of tracing, the sum of two terms 
Gra and Grr gives the full expression for conductance 
G{Ep), which have to cancel each other for \Ep\ > 6t 
in order to ensure vanishing of G{Ep) when tj^ — t is 
chosen. 

The result ior W — 7 disorder belongs to the non- 
perturbative transport regime. It is interesting, there- 
fore, to follow the behavior of Gra and Grr terms further 
throughout this regime, eventually reachiag Wc — 16.5 
where the whole band becomes localizedJl3 The generic 
LD transition point in three-dimensions (3D) is beyond 
NLSM treatment inasmuch as it corresponds to a strong 
coupling limit in the field theoretical language. Fur- 
thermore, recent descriptionCJ of Anderson localization 
in terms of an order parameter, obtained from a the- 
ory based on local approximation, suggest that probing 
the instability of the delocalized phase through WL-type 
corrections might be a daunting road to reach the LD 
transition in 3D. Figure ^ shows that around Wc the 
conductance G ^ 2e^/h is mostly defined by the Grr 
term — a situation completely opposite to the weak dis- 
order finding (W^ = 2 in Figs. l[and||). Therefore, at 
some intermediate disorder 7 <W < 10, within the non- 
perturbative transport regime, a transition occurs from 
Gra > Grr to Grr > Gra, ending up eventually with a 
case Gra < < Grr around the LD critical point. 

An inquisitive reader might have come up by now with 
the question as to what happens in the clean case. A 
mesoscopic ballistic sample attached,to two leads has 
non-zero point contact conductanceE3 of a purely geo- 
metrical origin since leads are widening into macroscopic 
reservoirs at infinity, where reflection occurs when the 
large number of conducting channels in the macroscopic 
reserjSapir matches the small number of channels in_the 
lead.Cj Such point contactpConductance is quaatizedjlj as 
a function of sample widthE3 or Fermi energyjiil which be- 
comes conspicuous when the number of quantized trans- 
verse propagating modes (i.e., the sample cross section) 
is small enough. A clean [sm = 0) toy sample 3 x 3 x 3 il- 



lustrates conductance quantization in Fig. ^, as obtained 
from the two-probe Kubo formula for G. In this case, 
both terms Gra and Grr are contributing in a non-trivial 
fashion to a stepwise conductance. 

What can be learned from these numbers is that 
Gra (|^) can serve as a decent approximation to the exact 
Kubo formula for the quantum conductance G only in 
very weakly disordered conductors. However, because of 
not conserving the current, the trace in Gra has to be 
performed over the whole disordered region (which is an 
enormous computational effort, and is therefore useless 
in the real-space computational practice). The essential 
outcome of this exercise is the explicit quantification of 
the difference between G and Gra- This, together with 
current conservation, are important points to bear in 
mind when using simplified expression Gra in analytical 
derwatipns and arguments of the quantum transport the- 
oryJ3ll3'E2l Finally, the strong disorder (nonperturbative) 
behavior of two different Green function terms, compris- 
ing the Kubo formula for quantum conductance, might 
provide a clue for the inadequacy of attempts to analyze 
genuine Anderson transition in 3D by probing instabil- 
ity of metallic phase to weak localization (perturbative) 
corrections. 
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